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We derive a classification and a measure of classical- and quantum-correlation of multipartite 
qubit, qutrit, and in general, n-level systems, in terms of SU(n) representations of density matrices. 
We compare the measure for the case of bipartite correlation with concurrence and the entropy of 
entanglement. The characterization of correlation is in terms of the number of nonzero singular 
values of the correlation matrix, but that of mixed state entanglement requires additional invariant 
parameters in the density matrix. For the bipartite qubit case, the condition for mixed state 
entanglement is written explicitly in terms of the invariant paramters in the density matrix. For 
identical particle systems we analyze the effects of exchange symmetry on classical and quantum 
correlation. 
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Quantum entanglement is an information resource; it 
plays an important role in many protocols for quantum- 
information processing, including quantum computation 
[l|, quantum cryptography 0], teleportation (3[, super- 
dense coding [i| , and quantum error correction protocols 
Q . Techniques for characterizing the bipartite entangle- 
ment and correlation of pure and mixed quantum states 
have enabled many advances in quantum information and 
the study of decoherence [|| . Many quantum information 
protocols use bipartite entanglement, but multipartite 
entanglement 0], also has quantum- information appli- 
cations, e.g., controlled secure direct communication [81 , 
quantum error correction Q , controlled teleportation [To| 
and secret sharing pd| . It has been shown that any insep- 
arable two-qubit states can be distilled to a singlet-state 
form with enough copies of the qubit-pairs [l2j and algo- 
rithms for multi-copy entanglement distillation for pairs 
of qubits have been developed (l3j . Moreover, multipar- 
tite entanglement offers a means of enhancing interfero- 
metric precision beyond the standard quantum limit and 
is therefore relevant to increasing the precision of atomic 
clocks by decreasing projection noise in spectroscopy fl4| . 
Here we use a representation of the density matrix for 
qubit, qutrit, and more generally, n-level systems con- 
taining 2, 3, . . . , and TV-parts, in terms of the correla- 
tions between the subsystems to quantify the classical 
and quantum correlation of multipartite systems. Our 
classification of correlation is in terms of the correlation 
matrix and its singular values, and our classification of 
entanglement of mixed states [l5( is associated with the 
Peres-Horodccki criterion fl9j , which we express in terms 
of additional invariant parameters in the density matrix. 
Separate measures of bipartite, tripartite, etc., correla- 
tion are required, since general mixed states can have 
bipartite correlation as well as higher subsystcm-number- 
corrclation. 



Werner [15( defined a mixed state of an TV-partite sys- 
tem as separable, i.e., classically- correlated, if it can be 



written as a convex sum, 

where p A is a valid density matrix of subsystem A, 
etc. Otherwise, Werner defined it to be entangled, i.e., 
quantum- correlated. Unfortunately, this definition of en- 
tanglement for mixed states is not constructive, since, in 
general, it cannot be used to decide whether a given den- 
sity matrix is separable or entangled. Moreover, a quan- 
titative measure of entanglement of multi-partite systems 
has proven to be difficult to devise. Note that studies of 
the best separable approximation to an arbitrary density 
have been carried out and have led to a proposal of a 
measure for entanglement (l6| . Furthermore, aspects of 
the geometry of separability and entanglement based on 
Schmidt decomposition have been studied and led to an 
analysis of the question of separability for the two-qubit 
case [13 ■ 

In what follows, we categorize classically-correlated 
and quantum-correlated states and characterize their cor- 
relation in terms the number of nonzero singular values 
[l8|], {di}, of the correlation matrix C, and characterize 
entanglement of bitpartite qubit systems using the Peres- 
Horodecki criterion [l9| which is reformulated totally in 
terms of the parameters used in forming the density ma- 
trix. 

First, let us consider a bipartite qubit system. For two 
uncorrelated qubits, call them A and B, we can write 
the density matrix as a product, pab = PaPb-, where 
the individual qubit density matrices can be written as 
pj = | (1 + nj • er /), where J = A,B, the crj are Pauli 
matrices for particle J and the Bloch vectors are nj = 
(er.j) = Trcr.jp j [20l| . For two correlated qubits, 

Pab = j [(1 + riA • cr A ) (1 + n B • cr B ) + cr A • C AB ■ cr B ] , 

(2) 

where the tensor C AB specifies the qubit correlations, 

G ij B = {<7i,A<Tj,B) - (<7i,A)(<7j,B) = (<7i,A<7j,B) -rH,AH>j,B- 

(3) 
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The density matrix pab is a 4x4 Hermitian matrix with 
trace unity, so 15 parameters are required to parame- 
terize it. The 3 components of n^, the 3 components 
of riB, and the 9 components Cy of the 3x3 matrix C, 
where we have no longer explicitly shown the subsystem 
superscripts, are sufficient for this purpose. 

Similarly for the bipartite qutrit case. The 3x3 den- 
sity matrix of a single qutrit can be written as p = 
I (l + |(Ai}Ai) where the Xi are the eight traceless Her- 
mitian Gellman matrices familiar from SU(3) (2lj . and 
(Xi) = Tr Xip. A bipartite qutrit density matrix can be 
parameterized in the form 

13 3 9 

PAB = ~ [(l+-(Ai,A)Aj,A) (l + -(Xj,B)X ] , B ) + -Xi 1 ACi j Xj, B \ 



Cij = (Xi t AXj t B) — (Xi,A)(Xj : B) 



(4) 
(5) 



where dj specifies the correlation between A^a and A^b. 
Here, pab is a 9x9 Hermitian matrix with trace unity, so 
80 parameters are required to parameterize it. The eight 
components of (X^a), eight components of (A^s), and 
64 components CV,- of the 8x8 matrix C are sufficient 
for this purpose. The same procedure can be used for 
bipartite 4-level systems using the 15 traceless 4x4 Her- 
mitian generator matrices for SU(4), and bipartite n- level 
systems with the n 2 — 1 traceless nxn Hermitian matri- 
ces. Likewise, a general qubit-qutrit 6x6 density matrix 
takes the form pab = | [(^+n i>A ^i,A) (1+f (Aj,b)A.,-,b) + 
j&i,ACijXj,B] with = (ct^aAj.b) - (cn,A)(Xj, B ), i = 
1, 2, 3 and j = 1, . . . , 8. 

Our bipartite correlation measure for an n-level and m- 
level system is based on the (n 2 — 1) x (m 2 — 1) correlation 
matrix C: 



4(n 2 < 



1) 



TrCC T 



4(< - 1) 



(6) 



states can be differentiated from classically-correlated 
states with 3 NSVs by applying the Peres-Horodecki 
(PH) partial transposition condition [l9| [which corre- 
sponds to changing the sign of n VtB and the matrix el- 
ements C^ B that multiply o~ Vt B in ([2]), and determining 
whether the resulting p is still a genuine density matrix 
— if it is, the state is classically correlated, i.e., unentan- 
gled but correlated] to the density matrices with 3 NSVs. 
The only categories that cannot be distinguished without 
use of the PH condition are the mixed-entangled and the 
classically correlated states with > 4 NSVs. 
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FIG. 1: Classification of two-qubit states. Categories can be 
experimentally distinguished by measuring xia, n_e, and using 
Bell measurements Q to determine the C matrix. 



Similarly, for a two qutrit pure state, the number of 
NSVs of C is zero for non-entangled states (the C matrix 
vanishes), three, if only two basis states are present in 
the entangled state, five, if one of the qutrits contains 
only two basis states but the other contains three, and 
eight if all three basis states are present. For classically 
correlated qutrit states, there are 1, 2, . . . , 8 NSVs for 
2,3,..., and 9 or more terms in the sum, etc. A similar 
classification in terms of the number of NSVs exists for 
qubit-qutrit and n-level systems. 

A general threc-qubit density matrix can be written as 



1 



'■./ 



PABC 



[(1 + n A ■ <ta) (1 + Hb • er B ) (1 + n c ■ er c ) 



r Tr (pab - PaPb) 2 



where n < = min(n,m). £ c = ^ _^ 
is a nonnegative real number. If C is a normal ma- 
trix TrCC T equals to the sum of the squares of 
its eigenvalues, but C need not be normal. £c is basis- 
independent; any rotation in Hilbert space leaves it un- 
changed. The normalization factor n^/^n 2 , — 1)] in ((6|) 
is such that the maximum possible value of £c is unity. 
£c measures both classical- and quantum-correlation. 
This measure of bipartite correlation was suggested in 
Ref. [22J for pure states and n — m. 

The correlation matrix C quantifies the correlation and 
the entanglement of bipartite states. For pure two-qubit 
states, the number of nonzero singular values (NSVs) of 
C is zero for non-entangled states (C vanishes), and three 
for entangled states. For classically-correlated states 
with two terms in the sum [see Eq. ([15])]. only one NSV 
occurs, two NSVs occur for three terms, three NSVs 
occur for four or more terms, and for entangled (i.e., 
quantum-correlated) mixed states there are three NSVs. 
These cases are summarized in Fig. [TJ Entangled mixed 



a a ■ C AB ■ cr B + er a ■ C AO ■ er c + (J B • C ao • rr c 



AC 



iBC 



/~] O' it A<Jj : B0'k : cDijk} 
ijk 



(7) 



where C AB , C AC , and C BC are the bipartite correla- 
tion matrices and the tensor that specifies the tripartite 
correlations is 

Dijk = (o-i,AVj,B&k,c) - (< J i,A)(o~j,B)(o'k,c) ■ (8) 
A tripartite qutrit state can be similarly parameterized: 

PABC = ^ [H(l + | J2(*i,l)Kl) + J J2 X i,lCij,U^,J 



27 



I, J i,j 



"jr \,AXj,BXk,cDijk,uK] 



(9) 



I,J,Ki,j,k 



3 



D ijk! u K = (KjXjjXk.K) - {Xi,i}(Xj,j)(X k ,K) • (10) 



Our tripartite correlation measure £d is based on the 
be written as 



correlation matrix D, Ed = K J^, j k -^'iik' which can also 



£d = KTr(p A Bc-PAPBPc- ^2 ^2 C if 1 Vi-JVj-j) 2 > 

(11) 

where Jv = 1/4 for qubits, and A' = 27/160 for qutrits 
with as replaced by As. £ B is also a basis- independent 
nonnegative real number; any rotation in Hilbcrt space 
leaves it unchanged. A tripartite system may have 
bipartite- as well as tripartite-correlation. The bipar- 
tite correlation of a tripartite system is the sum of the 
correlation for the three bipartite pairs, 



4(n 2 



]T TtC I " r (C I ' J ) T , (12) 
i,J{iitJ) 



where I, J = A, B, C. The density matrices of four- 
particle and higher qubit, qutrits, and n-level system 
states can be constructed similarly, but with increased 
complexity. For example, it is clear from Eq. how to 
generalize and obtain the four-particle correlation of four- 
particle systems: £e = A' ^ - . k l Efj k[ , where the four- 
particle-corrclation term of the four-qubit density matrix 
Pabcd is Y^ijkl a i,A° 'j,B<Tk,c&i,DEijki and A' = 1/8. 

We now present some examples of qubit and qutrit 
bipartite and tripartite correlated states. The maximally 
entangled bipartite qubit states are the Bell states, 
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on) ± nt)] 



i 

71 



[in> ± in)] • (13) 



For all these states, (<ta) = (o"b) = 0, i-e., = = 0. 
For the singlet, (ffi^ACj.B) = (the spins are oppo- 

sitely polarized). The density matrices of the Bell states 
are: 



Pm- 
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- cr A 
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+ er A 
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b) ■ 


(14) 



The correlation matrices of the Bell's states are diagonal 
and the correlation measure is Be = 1, i-e., they are 
maximally entangled. 

Let us now consider the Rashid pure states [23| . 

|0+) = (2cosh(20))- 1 /2 ( e -0| tt ) + e o U)) , (is) 



whose density matrix is p^ 



tanh(26»)cr Zi A] [Is - tanh(26»)o- ZiB ] + sech(20)(o- X)J 40- a;) . B - 
o- v ,ao- v ,b) + scch 2 (28) a z . A (Jz,b} ■ When 8 = 0, 



\4> + ) = | < f )+ ), and as 8 — >• ±00, an unentanglcd state re- 
sults. The nonvanishing correlation matrix elements are: 
C xx = sech(260,CV y = -sech(20), C zz = sech 2 (26»). 
Using ([5]) we obtain the correlation measure 
£c(|0 + )) = iTrCC 7, = !(2sech 2 (2(9) + sech 4 (26»)). 
The concurrence C [H HI is 

C(|0+» = V2(l-Tr [pa 2 ]) = scch(20) , 



PA = 



1 I 



tanh(20) 
1 + tanh(20) 



and the entanglement entropy is S = — Tr [pa log 2 pa]- 
These results are graphically presented in Fig. [2j All the 
measures equal unity for 8 = and decrease rapidly vs. 9. 




FIG. 2: (color online) Comparison of the Ec measure of cor- 
relation, the concurrence C and the entanglement entropy S 
for the Rashid pure states. 

Two-qubit classically-correlated states take the form 

P CC = ^^Pk {1 + n A ,k ■ <ta) {1 + n B . k ■ ct b ) , (16) 

fc>2 

with ^2 k Pk = 1 and pk > 0. The density matrix for the 
classically-correlated state can be written in the form of 
Eq. ([2]) with Bloch vectors 



n A = ^Pkn A ,k, ne = ^p k n B ,k 



and correlation matrix 
Cij = ^Pfc n lt A,k 

k 



n j,BM - / ]PI n i,B,, 



(17) 



(18) 



For example, for classically-correlated mixed states of the 
form p cc = (2sech 2 (20))- 1 /2( e -»|| t ) (| t | + e »| t ;) ^ 
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wc find that all the correlation coefficients vanish, except 
for C zz = — sech 2 (20), the density matrix in representa- 
tion (J2j is p cc = \ (l^ls - sech 2 (26») ct z> a(t z ,b), and the 
classical-correlation measure is £^P — | scch 4 (2#). 

It is elucidating to consider the Werner two-qubit den- 
sity matrix composed of a sum of a singlet state and the 
maximally mixed state, p w = p\^~) \ + ^j^l, or, the 
more general Werner two-qubit density matrix, 

p GW =p ^- ){q p- l + l^P 1 > (19) 



where = (2cosh(26»))- 1 / 2 (e - 

r educes to p w for 8 = 0. For p GW 
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iia = -ng = ptanh(20) z , 



(20) 



and 



^sech(26») 
sech(26>) 






1 - p + psech 2 (26l) 

m (21) 

The PH entanglement criterion [19( shows that this state 
is entangled if p[(l+2sech(26>)] > 1. Figure[3]plots the PH 
criterion limit and the correlation measure, £c{p,8) = 
J2i d 2 = 1 — p + (2p 2 + p)sech 2 (20), for the generalized 
Werner state. Note that the PH criterion is not obtain- 
able from C alone, but can be obtained using the invari- 
ant parameters £ = ^ <fj — n ^ A C n " B and • More 
explicitly, p[l + 2sech(2#)] = — £ + ^/£ 2 /4 - iia • n s , so 
the PH condition reads 



> 1 



(22) 



which can be written as the condition: the largest root of 
the quadratic equation, (x+£/2) 2 +£(x+£/2) + n.A-n-B 
ii. is greater than unity. Thus, mixed state entanglement 
is determined not only by C but by additional invariant 
characteristics of the density matrix, i.e., invariant char- 
acteristics composed of the parameters C, and 
used to form the density matrix (whereas the correlation 
is determined only in terms of C). The physical signif- 
icance of the scalar product 114 ■ ng as the projection 
of the expectation value of the spin of one qubit on the 
other, is clear, as is the physical significance of £ as a 
specific projection of the singular values of the correla- 
tion matrix that depends on the average spins and 
ng [25[. However, the physical significance of the PH 
entanglement criterion is not yet clear; i.e., the physical 
interpretation of Eq. (|22[) [or the quadratic equation] re- 
mains to be uncovered. But at least the PH condition is 
now expressed only in terms of the physical parameters 
appearing in the density matrix, rather than by the par- 
tial transposition condition, which is more removed from 
physical interpretation. 

As an example of a tripartite pure qutrit state, consider 
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FIG. 3: (color online) £c(p,0) versus p and 8 for the gener- 
alized Werner density matrix p GW , and the Peres-Horodecki 
entanglement criterion limit, p[l + 2sech(2#)] = 1, drawn on 
the p-9 plane and projected onto the Sc surface. 



where \v x ) = (1,0, 0) T , \v 2 ) = (0,1, 0) T , \v 3 ) = (0,0, 1) T . 
The qutrit bipartite correlation £q and tripartite corre- 
lation £0 are plotted in Fig. 01 The maximum of £d is at 
6\ = 8 2 = 0, where £jj = 1, but the bipartite correlation 
dips there. Three ridges of high bipartite and tripartite 
correlation occur, one at 8\ = 82 — negative, and two 
others at 120 degrees rotation from the first. 




(a) 




(23) 



FIG. 4: (color online) The bipartite and tripartite correlation 
measures, £c and £d for the tripartite-qutrit pure state (|23|) . 
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In most quantum information systems, qubits are dis- 
tinguishable, so there is no need to account for exchange 
symmetry, but for identical bosonic or fermionic systems, 
the density matrix must be properly symmetrized, e.g., 
Pab = $PabS- For two qubits, the symmetrization op- 
erator is S = |(1 + Pab) = § + a • &B and the anti- 
symmetrization operator is A — ^(1 — Pab) = j — \< J A • 
er B , hence, p s J^ = (| + \° 'A • cr B ) PAB (§ + \cr A • cr B ), 
Pab = (I ~ \ a A ■ &b) Pab (§ - j^a • <tb)- If the qubits 
are antisymmetric, pab = Pajs = -ApAB-A, the state 
must be pure singlet, p^g = j(l — a a • &b)\ it can- 
not be a mixed state, as opposed to p & J^ which can be 
mixed. If spatial degrees of freedom need to be included 
in the description, in addition to the internal degrees of 
freedom, a bipartite density matrix can always be writ- 
ten as a product of an internal (i.e., spin) part and an 
external (i.e., space) part. Hence, a symmetric density 
matrix p^g for the internal degrees of freedom must be 
multiplied by a symmetric [antisymmetric] spatial den- 
sity matrix for the spatial degrees of freedom {r^, r^} 



for bosons [fcrmions], and p A n g must be multiplied by 
an antisymmetric [symmetric] spatial density matrix for 
bosons [fermions] , so that the full density matrix has the 
right exchange symmetry. We show elsewhere that this 
has relevance to collisional shifts in atomic clocks [2(| . 

In summary, we have developed a classification of cor- 
relation for multipartite iV- level quantum systems by 
writing their density matrices in terms of SU(iV) gen- 
erators, and we defined a measure of correlation for such 
systems, based upon their correlation matrices. The en- 
tanglement involves not just the correlation matrix but 
also other invariant parameters in the density matrix for 
the system. This formulation can now be used for a vari- 
ety of applications, e.g., in the optimization of quantum 
gates and in the calculation of collisional clock shifts. 
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